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Abstract 

We present an extension of the two-point optical microrheology technique introduced by Crocker 
et al. [Phys. Rev. Lett. 85, 888 (2000)] to high frequencies. The correlated fluctuations of two 
probe spheres held by a pair of optical tweezers within a viscoelastic medium are determined using 
optical interferometry. A theoretical model is developed to yield the frequency-dependent one- 
and two-particle response functions from the correlated motion. We demonstrate the validity of 
this method by determining the one- and two-particle correlations in a semi-dilute solution of 
polystyrene in decalin. We find that the ratio of the one- and two-particle response functions is 
anomalous which we interpret as evidence for a slip boundary condition caused by depletion of 
polymer from the surface of the particle. 
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I. INTRODUCTION 



The dynamics of colloidal particles dispersed in viscoelastic polymer solutions are im- 
portant in both fundamental and applied science. Suspensions of particles in polymers are 
encountered in an vast range of chemical products including coatings, controlled-release 
drug formulations, personal care products, and filled polymer composites as well as being 
significant in many biological processes, such as intracellular transport inside the viscoelastic 
cytoplasm of a cell From a fundamental viewpoint, the characteristic lengths of polymers 
and colloids which differ by several orders of magnitude, suggest that these mixtures should 
be ideal candidates for coarse-grained descriptions. However while such approaches have 
been applied to equilibrium properties with considerable success there have been no compa- 
rable attempts, so far, to explain and quantify the rich dynamical and rheological properties 
of these complex mixtures. Understanding the dynamics of these complex systems remains 
an important goal of soft matter science. In recent years considerable experimental progress 
in the characterization of these systems has been made though the development of opti- 
cal microrheology techniques. The underlying idea, pioneered by Mason and Weitz 0, is 
that the Brownian fluctuations of a colloidal probe particle reflect the viscoelasticity of the 
medium in which the probe is embedded. By analysing the thermal motion it is possible 
to obtain quantitative information about the rheological properties of the polymer matrix 
over an extended range of frequencies not accessible to conventional rheometers. Optical 
microrheology experiments typically monitor the thermal fluctuations of a single probe par- 
ticle using a combination of video microscopy, optical interferometry [0, ^, or diffusing wave 
spectroscopy [0]. While microrheological techniques enjoy significant potential advantages 
over traditional rheological methods complications in interpretation have limited their ap- 
plication 1^. The problem is that the presence of a particle can perturb the surrounding 
medium (by, for instance, depleting polymer segments from near the particle surface) so 
that the dynamics of a single particle fails to reflect the macroscopic viscoelasticity of the 
medium P, 0. In effect, a particle diffuses in a "pocket" of material whose local rheological 
properties are not those of the bulk. To circumvent this limitation, Crocker et al. 0] con- 
jectured that the correlated fluctuations of two separated probe particles should measure 
the bulk rheological properties of inhomogeneous materials more accurately than traditional 
one-particle experiments. Using video microscopy they tracked the thermal diffusion of 
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pairs of probe particles embedded in a viscoelastic solution of guar and confirmed that the 
extracted viscoelastic moduli were in close agreement with conventional measurements. Sig- 
nificantly, the single particle results differed substantially, both in magnitude and scaling, 
from the two-particle results. Recently these observations have been confirmed theoretically 
by Levine and Lubensky |^. Using an analogy to classical electrostatics they have shown 
that while one-point techniques measure essentially a local rheology, two-point correlations 
determine the bulk rheology. Interestingly, these observations suggest that combined one- 
and two-point measurements could be a powerful tool to probe the nature of the specific 
interactions between particle and medium. 

In this paper, we present an extension of the two-point correlation technique to high 
frequencies using optical interferometry to track particles with high spatial and temporal 
resolution. Our approach allows access to a much faster range of timescales than previous 
measurements 0, as we show later. In consequence we are able to measure the rheological 
properties in a frequency range from typically ~ 10^ to lO^rad s~^. Accurate measurements 
of multi-point correlations are substantially more demanding than one-point measurements. 
To achieve adequate signal-to-noise levels we have been forced to use a more intense laser 
beam than usually employed in particle tracking experiments |^, We consider explic- 
itly the effects of optical-gradient forces and assume that each particle is held within a fixed 
harmonic potential. The two-particle correlations measured therefore differ from those dis- 
cussed by Crocker et al. [0 where diffusion occurs in the absence of a potential. We develop 
a generalised Langevin theory to account for the correlated time-dependent fluctuations in 
terms of the viscoelasticity of the medium. Our paper is organised as follows: In Section J|, 
we present an analysis of our two-particle microrheology experiments, and discuss how two- 
particle correlations are related to the viscoelasticity of the medium (Sec. [II B|) . Next, in 
Section |T| we describe our experimental setup and numerical procedures. Finally in Sec- 
tion 1^ we present measurements of one- and two-point fluctuations of two widely-spaced 
spheres immersed in a semi-dilute polymer solution which demonstrate the validity of our 
approach before we conclude in Section 0. 
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II. THEORY 



The fundamental concept underlying optical microrheology experiments is that the Brow- 
nian motion of a rigid spherical tracer particle is determined by the viscoelasticity of 
the surrounding medium. Typically the mean-squared displacement (MSD) (Ax^(t)) = 
{[x{t) — x(0)]^) of the particle trajectory is measured (here x{t) denotes the tracer particle 
position after a time t and the brackets denote a time or ensemble average). The MSD 
is related to the rheological properties of the medium though linear response theory. In 
general, the time- dependent frictional force F{t) on a sphere moving with a velocity u{t) 
through a viscoelastic medium may be written as 



where ^(t) is the instantaneous friction coefficient. The convolution on the rhs of Eq. |l| 
reflects the viscoelastic nature of the suspending fluid. Energy is stored as the Brownian 
particle diffuses, generating a "memory" of the particle's past motion. In consequence, the 
frictional force experienced at a time t depends on the particle velocity at all earlier times t'. 
This non-local time-dependence changes profoundly the Brownian trajectory so that from 
a measurement of the MSD the time- dependent friction coefficient C,{t) can be determined. 
The utility of microrheological measurements depend on extracting information about the 
rheological properties of the medium from the time dependence of C,{t). 

A. Correlated Brownian Motion in a Viscoelastic Fluid 

The optical gradient forces on a high refractive- index particle within a tightly-focused 
laser beam are well approximated by a harmonic interaction. Consequently, we model our 
experiments by analysing the dynamics of a pair of rigid Brownian particles of radius a, each 
held in a harmonic potential well, and separated by a distance r within a linear viscoelastic 
medium. Fluctuations in the positions of the two particles are correlated as a result of 
the hydrodynamic interactions which are transmitted through the viscoelastic matrix. The 
standard Langevin description (see for instance Refs. |1T|) for the Brownian motion of 
a pair of neutrally-buoyant particles of mass m is readily modified to include the effects of 
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the viscoelasticity of the medium, 
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Here the subscripted indices label the two Brownian spheres, k is the harmonic force 
constant of the optical traps (assumed identical), and //^(t) denotes the random Brownian 
forces acting on particle i {i = 1,2). The time-dependent frictional coefficients C,ij{t) in 
Eq. |] result from the generalisation of the viscoelastic "memory" (Eq. |l]) to the case of 
two interacting particles. The self term ^ii(t) details the force acting on one sphere when 
the same sphere is moving (and the second sphere is stationary) while the cross friction 
coefficient C,i2{t) describes the force generated on one sphere by the motion alone of the other 
sphere. Since the random thermal and frictional forces both have the same microscopic origin 
they must satisfy in equilibrium at a temperature T, the fluctuation-dissipation theorem 

{mt)f^{t')) = ksn^,it-t')m- 

The coupled generalized Langevin equation @ is simplifled by introducing collective 
and relative normal coordinates and using standard techniques to solve the resulting decou- 
pled equations of motion. Here we focus on the solution for the two-particle mean-squared 
displacement, namely 
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where the displacements are normalised by the equipartition result, (Axf) = 2kBT/k. The 
mutual two-particle mean-squared displacement dij{t) for i j details the correlation be- 
tween the thermal fluctuations of the two separated particles. The self term dii{t), is the 
MSD of particle i in the presence of a second particle. The exact solution for the two-particle 
displacement function is from Eq. 0, 
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where the tilde denotes a Laplace transform, f{s) = f{t)e~^^dt, s is the Laplace variable 
and ^± = ± ^12- The complete solution (H) is still too complicated to use to analyse 
experimental data. In order to derive a simpler result we make two simplifying assumptions. 

We start by considering the importance of particle inertia. The mass m of the particle 
determines the decay time tb ~ m/GTrrjocCi of the particle velocity autocorrelation function. 



Here ?7oo is the high frequency viscosity of the medium. For frequencies small compared to 
ujb = 1/tb the momentum of the Brownian particle will have relaxed to zero and we may 
neglect particle inertia. In our experiments, uob is of order 10 MHz so neglecting particle 
inertia is an excellent approximation at frequencies of experimental interest {uj < 20 kHz). 

Secondly, we consider only the situation where the two probe spheres are separated by a 
distance r which is large in comparison to the sphere radius a. In this limit, Batchelor 



has argued from general symmetry considerations that the hydrodynamic functions depend 
on the dimensionless sphere separation p = r/a as, 

= (5) 

P 

where the coefficients and bi are constants for a given medium. This asymptotic solution 
reveals that provided the two spheres are far enough apart, the self friction ^.u is, to leading 
order, independent of the separation p. The mutual friction coefficient ^12 by contrast 
depends sensitively on the pair separation varying as 1/p at large r so that the ratio ^2/^11 
is a small parameter which shrinks with increasing separation as 1/p. Consequently the 
exact expression for the two-particle displacements (Eq. ^ may be simplified by expanding 
dij in powers of ^12/^11- At large p, only the constant terms and the part which decays with 
distance as 1/p contribute to the displacement field. Retaining these leading terms yields 
the desired simplified expression for the correlated sphere displacements. Ignoring particle 
inertia, the one- and two-particle displacements may then be expressed as 

k 
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where i ^ j- Rearranging these equations gives expressions for the friction coefficients in 
terms of experimental quantities, namely 
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Measurements of the single-particle fluctuations yield the self friction ^ii{s), whereas the 
mutual friction coefficient is found from a knowledge of both single-particle and two-particle 
positional fluctuations. Our analysis is valid for any viscoelastic media, we make no as- 
sumptions about the homogeneity of the medium. We assume only that: (i) the particle 
pair separation r is large in comparison to the sphere radius a, and (ii) particle inertia is 
unimportant. 



B. Relation to Rheology 

The relationship between the friction coefficient C,ij{s) and the viscoelasticity of the 
medium is a challenging problem in fluid mechanics. Some of the complications are il- 
lustrated in Figure which shows schematically the microstructure of a colloid-polymer 
mixture. If the polymer is non-adsorbing then around each particle is a depletion zone of 
width ( where the polymer concentration is less than in the bulk. Physically, we expect the 
motion of an individual sphere to depend on the local environment rather than the rheology 
of the bulk solution. So it seems reasonable to expect the depletion zone around a particle 
to speed up diffusion, in comparison to the case where the medium is homogeneous. 

While such specific interactions between probe and medium have not been considered, 
Levine and Lubensky have analysed in detail the simpler case of a viscoelastic continuum. 
Over a wide range of frequencies, they predict that the single particle friction coefficient 
satisfies the generalised Stokes-Einstein relation, 

~ 67iaG{s) 

^ii{s) = , (8) 

where G{s) is the Laplace shear modulus of the homogeneous medium. This result reduces 
to the well-known Stokes expression = Gnrja) for the friction of a rigid sphere in a viscous 
fluid IP] if the Laplace modulus is replaced by its Newtonian limit, G{s) = srj. Deviations 
from Eq. ^ are expected if the local environment around each sphere is not the same as the 
bulk. 

While probe-particle interactions have a signiflcant effect on single particle motion, 
Crocker et al. have suggested that they should have a much weaker influence on the 
correlated fluctuations of pairs of particles. This conjecture has recently been confirmed 
theoretically |H]. In a homogeneous viscoelastic medium, the two-particle friction coefficient 
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has the form, 

ps 

This equation is vahd provided the pair separation r is large in comparison to the sphere 
radius a. In the Newtonian hmit, we recover the asymptotic form of the mutual friction 
coefficient for a viscous fluid (^^ = —dnrja/p [0). 

III. MATERIALS AND METHODS 
A. Dual-beam Optical Tweezers 

The details of the dual-beam optical tweezer equipment used in the present experiments 
have been described in a preceding publication |11[]. Figure ^ shows schematically the 



apparatus. In brief, two spherical colloidal particles are held by optical gradient forces near 
the focus of a pair of orthogonally-polarised laser beams (A = 1064 nm). The mean separation 
r between the probe particles was varied by altering the positions of the trapping lasers with 
an external computer-controlled mirror. The thermal fluctuations in the position of each 
probe particle were monitored by recording the interference between the forward scattered 
and transmitted infra-red beams with a pair of quadrant photodetectors. Custom-built 
current-to-voltage converters allowed the in-plane positions of both spheres to be recorded 
with a spatial resolution of ~ 1 nm at time intervals of 50 /is. 

The optical gradient force on a sphere displaced from the focus of the laser beam is 
harmonic. The force constant k of each trap was determined in a separate calibration 



experiment Colloidal spheres dispersed in decalin were trapped in each of the two 

beams in turn. The single particle mean squared displacement (Ax^(t)) was measured and 
fltted to (Ax^(t)) = (Ax^) [1 — exp(— t/r)] to yield the trap stiffness k = GTcqa/r from known 
values of the particle radius a and solvent viscosity rj. The intensity of the two beams was 
adjusted until the stiffness of the traps differed by less than ~ 5%. For a given geometry 
and laser intensity, the harmonic force constant k is a function only of the refractive index 
mismatch between particle and medium, so these calibration constants are expected to hold 
for the polymer data presented below. The root-mean-square fluctuation in the particle 
position in each trap was ca. 50 nm. 
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B. Numerical Methods 



The normalised one- and two-particle mean-squared displacements dij (t) were calculated 
from the measured particle trajectories {xi{t),X2{t)} using fast Fourier transform algorithms 
Tl| . The Laplace transform of the one-particle MSD dii{t) was determined from a regularised 



fit to a linear superposition of exponential terms, 



N 



du{t) = l-Y,Lke-"^\ (10) 

k=l 

where N is the number of terms, is the decay time specified on a logarithmic grid and the 
coefficients Lk were determined by requiring that the spectrum of relaxation times remain 



smooth and penalizing sums proportional to their second derivative |T5[. The curvature 
penalty suppresses unphysical wild oscillations in the solution. The constrained regularised 
solution to Eq. |T0| was found using the CONTIN algorithm [|16] . Once have been deter- 
mined the Laplace transform may be calculated exactly as, 

N 



di 



--Vl,^^. (11) 



The Laplace transform of the two-particle MSD dij(t) was evaluated by a two step nu- 
merical procedure. First, dij{t) was fitted to a series of cubic polynomial splines. The 
contribution to the transform from each time interval was then expressed in terms of a se- 
ries of incomplete gamma functions which were evaluated numerically. Tests showed that 
errors in the transform were of the order of 5%, except near the frequency extremes where 
truncation errors became more significant. 

IV. RESULTS 

To test our approach we have made measurements on a viscoelastic semi-dilute solution of 
polystyrene (PS) with an average molecular weight of 10^ (Mw/Mn = 1.19) in a mixed cis- 
and trans-decalin solvent (volume fraction czs-decalin 0.48). Linear polystyrene in decalin 
is a well studied system 1^^. The B-temperature of PS in an equal volume mixture of cis- 
and trans-decalin is ~ 16 °C so at the temperature of our measurements (~ 23 °C) decalin 



remains a near-B solvent. Using literature values |17| we estimate the polymer radius of 



gyration at 23 °C as Rg ~ 102 nm, slightly larger than the value under B-conditions of 87 nm. 



The overlap concentration c* is 3.7 mg cm~^. Sterically-stabilised poly(methyl methacrylate) 
spheres of radius a = 0.643 /im were added as probe particles at a volume fraction (p ^ 10^^. 
The viscous suspensions were loaded into flat, rectangular capillary tubes 170 /im thick which 
were hermetically sealed with a fast-setting epoxy glue and mounted onto a microscope slide. 
The samples were left for at least forty minutes to allow the particles to settle to the bottom 
surface before the beginning of an experiment. The trajectories of two spheres were collected 
at ten roughly even-spaced pair separations between 3 /im and 17 fim. 2^^ data points were 
collected at each separation at a sampling frequency of 20 kHz. To avoid wall effects, we 
analysed only spheres at least 15 /im away from the capillary walls. All measurements were 
performed at room temperature. 

The time-dependent trajectory of a sphere diffusing in a PS solution is illustrated in 
Fig. ^. In contrast to cross-linked gels, the entanglements present in a semi-dilute polymer 
solution are only temporary so on a sufficiently long time scale the solution is expected to 
be purely viscous. Consequently the trajectory, when averaged over a long enough period, 
will depend only on the external potential and not on the viscoelasticity of the solution. 
The optical forces are harmonic so that at long times the particle probability distribution 
will be Gaussian and the mean squared displacement will plateau at the equipartition limit, 
(Ax^(t)) = 2kBT/k. This behaviour is seen in Fig. |^ (d). At short times, the particle statis- 
tics are also Gaussian although diffusion is now hindered by the temporary entanglements 
found in the surrounding semi-dilute polymer solution (Fig. ^ (a)). The non-Gaussian dy- 
namics at intermediate times (Fig. |^ (b) and (c)) are most intriguing. The trajectory is 
highly anisotropic and does not resemble the isotropic diffusion seen on either shorter or 
longer time intervals. The particle seems to make infrequent large "jumps" over a distance 
which is comparable to the polymer radius of gyration. This behaviour suggests that the 
instantaneous micro-environment around each sphere in a polymer solution is heterogeneous 
and evolves in time. 

To characterise the mechanical environment in detail we calculated the mean-squared x- 
displacement , (Ax^(t)), from the trajectory of a single particle. The data was recorded from 
a pair of widely-separated spheres trapped within a PS solution. Fig. ^ depicts one of the two 
MSDs, normalised by the plateau displacement, at a number of different pair separations. 
Data for the second sphere was essentially identical and is not shown. As expected, there 
is no systematic dependence of the single particle MSD on the pair separation r, provided 
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the two spheres remain well separated (r ^ a). The measured (Ax^(t)) is sensitive only 
to the local rheological environment around each particle. Because of the high frequency 
elasticity of the polymer solution, {Ax'^{t)) increases approximately as t^'^ at early times 
rather than linearly, as expected for diffusion in a purely viscous medium. At long times 
(Ax^(t)) reaches the equipartition plateau. 

Using the same trajectory data, we have also calculated the two-point correlations. We 
expect that for widely-spaced particles the strength of the correlated positional fluctuations 
will decay inversely with the particle spacing. To confirm this, we plot in Fig. ^ the mutual 
mean-squared displacement dij(t) multiplied by the dimensionless spacing p. The data 
measured for a range of different pair spacings is seen to collapse onto a common curve over 
a wide range of times, although deviations are evident at very long times. The variations 
seen at long times do not depend systematically on the pair spacing p and are probably a 
reflection of poor statistics or mechanical vibrations. The peak seen in dij (t) is a consequence 
of the applied optical potential since for very long times the correlations between the two 
particles must decay ultimately back to zero. 

The one- and two-particle mean-squared displacements contain features which originate 
from the viscoelasticity of the polymer solution and from the applied optical forces. To 
separate the two contributions we use the analysis of Sec. ^ Inverting duit) and dij(t) 
yields friction coefficients which, apart from a simple scaling, do not depend on the applied 
potential. The computed friction coefficients are shown in the inset in Fig. ^. The frequency 
dependence seen is in sharp contrast to that expected for a purely viscous fluid where 

remains independent of frequency. The steady decrease in the particle friction with 
increasing frequency reflects the viscoelasticity of the polymer solution and in particular 
its ability to store energy elastically at high frequencies. To discuss further the connection 
between the particle friction and the medium rheology we focus on the frequency scaled 
functions sC,ii and —sp^ij. In the continuum model of Levine and Lubensky |l^ both 
functions are predicted to be simple multiples of the Laplace-transformed shear modulus 
G{s) so that the ratio —sp^ij/s^a should be frequency-independent and from, Eqs. H and 
^, equal to 3/2. The experimentally-derived values for s^a and —sp^ij are shown in Fig. |^. 
In partial accord with these predictions, we find that the two scaled friction coefficients do 
exhibit a very similar functional form (except possibly at the highest frequencies). However 
the ratio —sp^ij/s^a is not 3/2 but of order unity. 
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While we do not have a complete understanding of this discrepancy one possible explana- 
tion is that the nature of the shear coupling between the particle and the medium differs in 
the one- and two-particle situations. In the immediate vicinity of a probe particle a depletion 
zone exists with a reduced polymer segment concentration. While hydrodynamic forces will 
still couple the motion of the probe particle to the polymer matrix, the depletion zone could 
modify the nature of the boundary conditions at the particle surface. For a rigid spherical 
particle, the choice of boundary conditions changes the friction coefficient by a factor of 3/2 



from Q-Kfja for stick to Airrja for slip 0. Cardinaux et al. [18| have recently suggested that 
a slip, as opposed to a stick, boundary condition might apply to single particle motion in 
giant-micellar solutions. Assuming the depletion of polymer has a similar effect, then we 
should identify ^^(s) with 47raG'(s)/s rather than Eq. |]. In the case of slip, the two-particle 
friction coefficient ^ij is —p^ij{s) = 4:7raG{s)/s fl^. In this case the ratio —p^ij/^u is 1.0 
which is in close agreement with the value found experimentally. Further experimental and 
theoretical work is necessary to confirm this picture. A detailed comparison between the one- 
and two-particle friction coefficients as a function of polymer concentration and molecular 
weight is currently underway and will be the subject of future publications. 

The frequency-dependent complex shear modulus G*{uj) = G'{uj) + iG"{uj) may be de- 
termined from either the one- (after allowing for the slip boundary conditions) or two-point 
friction coefficients. Fig. |^ shows the storage and loss moduli computed from tweezer mea- 
surements for a 2.9 c* solution of polystyrene. The data exhibits the features expected for 
an entangled flexible polymer solution [|lOl. At high frequencies, above the maximum re- 
laxation time of the solution, the shear and loss moduli are expected to show a power law 
dependence on uj. Under B-conditions, G' and G" are predicted by the Zimm model to vary 
as uj"^/'^. Fig. 1^ shows that at high frequencies the complex modulus G*{uj) follows closely the 
predictions of Zimm theory. At low frequencies the elasticity of the entanglement network 
contributes and the measured G'{uj) is increased above the Zimm predictions as expected 
for a semi-dilute solution. 



V. SUMMARY 

In this paper we have demonstrated that combined measurements of one- and two-particle 
thermal fluctuations can be used to provide new insights into the the nature of the interface 
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around a probe sphere embedded in a viscoelastic medium. We have measured the correlated 
Brownian motion of two widely-separated spheres held by a pair of tightly-focused laser 
beams within a polymer solution. Photodetection of scattered light was used to record the 
trajectories of the trapped particles with high spatial and temporal resolution. This new 
technique differs from previous two-point microrheological studies [0] in its spatial resolution 
and access to high frequencies; frequencies up to lO'^rad s~^ are readily detectable. In 
comparison, the frequency range of previous two-point measurements [0] have been limited 
by the speed of video microscopy to ~ lO^rad s~^. Faster A/D data acquisition should 
enable the optical tweezer technique to be extended to still higher frequencies approaching 
lO^rad s-\ 

In order to achieve good statistics on two-particle correlations we have used a significantly 
higher-powered laser beam than previous (single-particle) optical tracking studies |jl|, 0, 
As a result we are no longer able to ignore the effects of optical forces applied by the laser 
beam on the measured particle fluctuations, as previous studies have assumed implicitly 

1^, ||. We account for the optical-gradient forces by analysing the dynamics of two 
harmonically-bound particles embedded in a general viscoelastic medium. By making two 
simplifying assumptions, which we show are valid in our experiments, we develop simple 
expressions for the frequency- dependent single- and two-particle friction coefficients, in terms 
of measurable quantities. We expect the single-particle friction coefficient to be sensitive 
to the rheological properties of the interfacial zone immediately surrounding a probe sphere 
while two-point measurements should allow the bulk rheological properties to be determined. 

We have demonstrated the validity of our approach by analysing the one- and two-point 
correlations measured for 1.28 fim diameter poly (methyl methacrylate) spheres suspended in 
a semi-dilute solution of polystrene. One- and two-particle friction coefficients are computed. 
We confirm that over a wide range of particle separations, the two-point correlations vary in- 
versely with the particle spacing r while the one-point correlations are independent of r. The 
two friction coefficients show a very similar dependence on frequency with the single-particle 
friction coefficient ^m(s) approximately equal to the corresponding two-particle friction co- 
efficient —p^ij{s). While viscoelastic theory confirms that the two functions should have a 
similar frequency dependence, the ratio measured is in quantitative disagreement with recent 
theoretical calculations for a viscoelastic continuum which predict a ratio of 3/2. We 
suggest that this discrepancy reflects the existence of a depletion zone around each particle 
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in a non-adsorbing polymer solution. As a consequence, the shear coupling between particle 
and medium may be approximated by a slip rather than a stick boundary conditions on the 
particle surface. With this assumption, either the one- or two-particle friction coefficients 
can be used to obtain the bulk rheology. 
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FIG. 1: Schematic illustration of two-point microrheology. Two spheres of radius a are suspended 
a distance r apart in a semi-dilute polymer solution. The non-adsorbing polymer is depleted from 
an interfacial region of width surrounding each particle (shown hatched). The dynamics of an 
individual particle is a sensitive function of the width and properties of the interfacial zone. By 
contrast, the correlated pair fluctuations depend upon the mechanical properties of the medium, 
averaged on the longer scale r. In the limit when r 3> a 3> C the pair fluctuations measure the 
bulk rheology. 
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LI L2 

FIG. 2: Dual-beam optical tweezer setup: Optical traps are generated by focusing two orthogonally- 
polarized IR laser beams (LI and L2) using a high numerical aperture microscope objective. Inten- 
sity shifts caused by interference between the direct beam and light scattered by the trapped sphere 
are imaged onto a pair of quadrant detectors (QDl and QD2). A polarizing beam splitter (PBS) 
is used to separate the two orthogonal signals. The x and y positions are derived by combining the 
voltages from the four segments of the quadrant photodiodes using low-noise analog electronics. 
The signals are then digitized and stored on a PC. 
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FIG. 3: The measured two-dimensional trajectory of one of two widely separated spheres in a 
semi-dilute PS solution (c = 1.7 c*) over a duration of (a) 0.0128 s, (b) 0.0512 s, (c) 0.2048 s, and 
(d) 0.8192 s. The sampling frequency is 20 kHz. 
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FIG. 4: Single particle mean squared displacements measured for a pair of widely-separated par- 
ticles (c = 1.7c*). The different curves depict data collected at three different pair separations 
ranging from p = 4.53 to 16.96. No systematic variation with p is evident. 
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FIG. 5: The time dependence of the scaled mutual two-particle mean squared displacement pdij{t) 
measured at four different pair separations p (c = 1.7 c*). For intermediate times all of the data 
collapse onto a common curve showing that the strength of the two particle correlations varies 
inversely with the pair spacing p. The solid line is the average for all of the data collected. 
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FIG. 6: Comparison of the frequency dependence of the scaled two-particle —{s/k)p^ij{s) (solid) 
and single-particle friction coefficients {s/k)^ii{s) (dashed) determined from the correlated fluc- 
tuations of two widely-spaced particles in a semi-dilute (c = 1.7 c*) polymer solution. The inset 
diagram shows —{p/k)^ij{s) (solid) and S,n{s)/k (dashed) as a function of frequency. 
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FIG. 7: Scaled storage QiraG' {uj) jk (open circles) and loss moduli 6TraG"{u!)/k (filled circles) as 
function of frequency for 2.9 c* solution of polystyrene in decalin. The dashed lines show predictions 
from Zimm theory. 
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